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Abstract. The Lockman Hole is well known as the region with the lowest neutral atomic hydrogen column density 
on the entire sky. We present an analysis of the soft X-ray background radiation towards the Lockman Hole using 
ROSAT all-sky survey data. This data is correlated with the Leiden/Dwingeloo survey (Galactic Hi 21cm-hne 
emission) in order to model the soft X-ray background by using radiative transfer calculations for four ROSAT 
energy bands simultaneously. It turns out, that an important gas fraction, ranging between 20-50%, of the X-ray 
absorbing material is not entirely traced by the H I but is in the form of ionized hydrogen. Far-ultraviolet absorption 
hne measurements by FUSE are consistent with this finding and support an ionized hydrogen component towards 
the Lockman Hole. 
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1. Introduction 

In this paper we focus on the Galactic plasma emission 
with emphasis on the coronal gas in the Milky Way halo. 
We concentrate our study on the Lockman Hole region 
on the no rthern high Galactic la t itude sky. The Loc kman 
Hole (e.g. iLockman et"al] l)l986|) . Ijahoda et all (TogO)) is 
the sky region with the absolute lowest H i column density 
(iVni ~ 4.5 • 10^^ cm^^). Therefore, it is considered to 
be the most transparent window to the distant Universe. 
To test this assumption, we investigate the photoelectric 
absorption of X-rays by the interstellar medium (ISM). 

Towards this aim, we performed a correlation anal- 
ysis of the ROSAT all-sky survey data with the 
Leiden/Dwingeloo Hi 21 cm-line survey. Both surveys are 
the state-of-the-art data bases for such an approach. In 
the soft X-ray energy regime (E < 1 keV) the mean free 
path length of the X-ray photons are only a few tens of 
parsec at maximum, assuming a density of n ~ 1 cm~^ 
ijMcCammon fc SandersI [l99o|) . The exact knowledge of 
the amount of photoelectric absorbing material is essen- 
tial to evaluate the soft X-ray radiative transfer through 
the Galact ic ISM. Previous inves tigations of the Lockman 
Hole area ijSnowden et alJir994a|) showed that diffuse soft 
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X-ray emission originates beyond some Hi gas clouds 
on the line of sight. They concluded that in general the 
1/4 keV X-ray emission appears to be anti-correlated with 
the Nui. In detail, significant deviations between the lo- 
calization of individual H i clouds and the X-ray intensity 
minima are suggested. Moreover, ISnowden et al.l l)l994a|) 
proposed that the X-ray emitting halo g as located beyond 
the Draco cloud fenowden et aljri99l|) and the one ob- 
served towards the Lockman Hole are different in temper- 
ature T and emission measure indicating a different ori- 
gin of both plasmas. This supports the picture of a patchy 
Milky Way halo gas. 

In this paper, we analyze a much larger portion - en- 
closing the Lockman Hole - of the northern high Galactic 
latitude sky. We use the ROSAT all-sky survey data of four 
energy bands, sampling the soft energy regime: 0.25 keV 
< E < 0.82 keV. The Leiden/Dwingeloo Hi survey data 
serves as a tracer for the photoelectric absorption, eval- 
uated by the H i column density along the lines of sight. 
We utilize the velocity information inherent to the H i line 
spectra to unravel spatially the composition of the ISM 
towards the Lockman Hole. In particular, we focus on 
the role of the Warm Neutral Medium (WNM) and the 
intermediate-velocity cloud gas (IVC) towards this posi- 
tion on the sky. 
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The paper is organized as follows: Section 2 presents 
the X-ray and Hi data. Section 3 outlines our approach 
to correlate both data sets, summarizes the best fit pa- 
rameters for the temperatures, intensities and absorbing 
column densities. Section 4 describes the far-ultraviolet 
absorption line analysis of the Lockman Hole. In Section 
5 we present our conclusions. 

2. The Data 

2.1. Hi Data 

The Le iden/Dwingeloo survey of Galactic neutral hy- 
drogen l|Hartmann &: BurtorJll997(l covers the entire sky 
north of (5 = —30° with a true-angle spacing of 0?5 in 
both I and b. The spectrometer used has 1024 channels 
and a velocity resolution of 1.03 km s^^ which leads to 
a LSR velocity coverage of -450kms"i to +400kms~^ 
Hence, most high-velocity cloud (HVC) emission is encom- 
passed by this velocity interval. The rms-noise of the mea- 
sured brightness-temperature amounts to ATb = 0.07 K. 
Moreover, the spe ctra are corrected for the influence 
of str ay radiation l|Kalberla et al.lll980t iHartmann et alJ 
Il996|) . which is of greatest concern in galactic Hi-surveys. 

We extracted the sky region bordered by 6 = 12°. . . 75° 
and I — 99°. . . 166° including the Lockman Hole. The 
lowest Hi column density region ijLockman et alJ Il986l 
(A^Hi ~ 4.5-10i3cm-2)) is focalized at {l,b) = (150°, 53°). 
According to the large stretch in Galactic latitude, the 
contrast between the highest column density region at low 
Galactic latitude and the Lockman Hole is about a factor 
of 25 (equivalent ~ 14 dB). This high column density con- 
trast gives us the opportunity to disentangle the contribu- 
tion of the different diffuse X-ray emission components, by 
focussing on regions where each individual X-ray emission 
component is observable with high significance. The H i 
data are binned into velocity integrated maps which are 
projected on the ROSAT grid for a pixel-by-pixel correla- 
tion. This is described in Sect. 13.31 The pixel size is set to 
48', about twice the ROSAT all-sky survey data analyzed 
here and about \/2 of the Leiden/Dwingeloo beam size. 



2.2. X-ray Data 

The X-ray data of the Lockman are a were extra cted from 
the R OSAT all-sky survey (RASS) (jSnowden et al. 1995t 
Il997|) . The data used have an angular resolution of 24'. 
We investigate four ROSAT energy bands: the Rl-band 
{E ~ 0.25 keV), the R2-band {E ~ 0.31 keV), the C- 
band (E - 0.28 keV) and the M-band (E ~ 0.82 keV) 
l)Snowden et alJll994tJ). The data are corrected for scat- 
tered solar X-ravs llSnowden fc Frevbere 19931. particle 
background ((Plucinsk y ^1^ ^l.| ll993ll. as well as long-term 
X-ray enhancements IjSnowden et alJll995|l . X-ray point 
sources wer e removed down to a minimum count rate of 
0.02 ctss"^ (ISnowden et al.lll997|) . The statistical signif- 
icance of the X-ray data were evaluated by the ROSAT 
uncertainty maps which only account for the number of 



photon events: they do not include any systematic uncer- 
tainties produced by any non-cosmic X-ray events. 

3. Method 

3.1. Radiative transfer equation for the ISM 
3.1.1. The effective photoelectric absorption 

In the soft X-ray energy regime, the photoelectric ab- 
sorption is the dominant process of X-ray attenuation. 
According to Morrison & McCammon (1983), the photo- 
electric absorption cross section is a strong function of 
energy {a oc E~^) and amounts to some a ~ 10"^" cm^ 
(according to ROSAT C-band) for an ISM of solar metalic- 
ity and n ~ 1 cni^"^. Both findings highlight two important 
facts: a) the lower the X-ray photon energy the higher the 
probability of photoelectric absorption and b) all lines of 
sight through the Galactic ISM are optically thick for soft 
X-ray emission - except for the Lockman Hole. As pointed 
out in Sect. 12.21 we analyze the X-ray intensity distribu- 
tion of broad energy bands. According to this approach, 
it is necessary to evaluate the X-ray band-averaged pho- 
toelectric absorption cross section, the so-called effective 
photoelectric absorption cross section: 



J_ j JSiE)-TiE)-e-'^'(^y^-'-dE \ 
7Vh/M J S{E)-e—'(Ey^^'-dE 



Here S{E) denotes the spectrum of the X-ray source, T{E) 
the telescope and detect or response function, cr'{E) the 
iMorrison fc McCammonI l(l983') photoelectric cross sec- 
tion, A^tt the obs erved H i column density (for details see 
iKerp et al.l ll99jj^. Using this effective photoelectric ab- 
sorption cross section (t(A^hi) yields an inherent ambigu- 
ity: without an educated guess of S{E) it is not possible 
to determine cr(A^Hi). In consequence, we have to compare 
the observational data with a variety of different spectral 
models to overcome this ambiguity. 

3.1.2. Galactic iYni 

The distribution of the general ISM of the Milky Way 
is traced by the Galactic Hi 21 cm- line emission. It is 
thought that H i - as the most abundant element in space 
- traces the distribution of the metals quantitatively. 

Early Hi surveys already revealed that Hi clouds 
exist next to the rotating Galactic disk, which can- 
not be modeled within the framework of co-rotating 
ISM. Based on these early findings, today we distin- 
guish between low-, interrn ediate- and high- velocity clouds 
i Dickev fc Lockrria,nlll99fl|l . In the following we will show, 
that low- and intermediate-velocity clouds can be iden- 
tified as soft X-ray shadows ag ainst the distaii t diffu se 
X-ray background emission (e.g. 'S nowden et alJ Emh ^. 
To differentiate between these different constituents of 
the ISM we bin the data into velocity integrated maps 
covering IwlsrI < 25kms~^, |wlsr| < 75kms~^ and 
I^^lsrI < 100 kms~-^. We show that the |wlsr| < 75kms~-^ 
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already covers Hi gas belonging to the M81/M82 group 
of galaxies. This intra-group Hi gas allows a critical test 
of our X-ray vs. Hi correlation analysis because it is acci- 
dentally included in the radiative transfer equation as an 
absorber within the Ii^lsrI < 75kms~-^ map of the Milky 
Way halo. 

3.2. The individual X-ray source terms 
3.2.1. The Local Hot Bubble 

Early sounding rocket experiments as well as the 
Wisconsin X-ray sky survey established the existence of 
hot X-ray emitting plasma within the local environment of 
the Sun. This c oronal gas fills the i rregularly shaped local 
void of matter ijMcCammon fc San ders..,199Q:l - frequently 
called the Local Hot Bubble (LHB). Therefore, the X-ray 
intensity of the Local Hot Bubble varies appreciable across 
the entire sky: /t.hb = (2.5 . . . 8.2) • 10~^ cts s~^ arcmin"^ 
I Snowden et alJ Il998t) . /lhb is the first term in our ra- 



diative tra nsport eq uation (see Eqn.[21). Considering the 
results of .Sfeir et alJ (1999) wc know that only a few local 
clouds are intermixed with the hot X-ray emitti ng plasma 
of the Local Hot Bubble fe.g. iKerp et al.lfl993|) . Because 
of this finding, we assume that /lhb is not absorbed by 
any intervening gas along the line of sight. 



3.2.2. The Milky Way halo 

The existence of hot coronal gas within the Milky Way 
halo was one of the major discoveries of the ROSAT 
mission. Moreover, it was found that the X-ray emis- 
sion from the LHB has a superposed diffuse comp onent 
llSnowden et alJll99lHHerbstmeier et al..,1995: Kcrp et alJ 
Il99fil IPietz et al.lll998l) in excess to the contribution of 
the extragalactic background radiation. It is beyond any 
doubt that diffuse X-ray emission originates within the 
Milky Way halo, but it is a matter of debate as to how 
many thermal pla sma component s arc localized within 
the Galactic halo l|Kuntz fc Snowden .2000') . Their corre- 
sponding individual temperatures and emission measures 
are also a matter of debate. Here we follow the approach 
of a hot single halo component which will turn out to 
be sufficient. We vary the temperature of the Galactic 
halo plasma and search for a best fit temperature simul- 
taneously consistent with all four ROSAT energy band 
data. The Milky Way halo is represented as /halo in the 
radiative transfer equation (see Eqn.|21). This diffuse X- 
ray emission component is absorbed by the Galactic ISM 
along the line of sight. Following the latest results on the 
scale height of the highly ionized species l| Savage et al.l 
Il997. 2003) and prior Hi vs. X-ray correlation analysis 
l|Pietz et al.lll99al) . the s cale height of the Milky Way halo 
plasma is about 4kpc fKalb erla fc Kerdll998l) . We can 
safely assume that all ISM components with a lower scale 
height, the Warm Neutral Medium (WNM), the Warm 
Ionized Medium (WIM) and the Cold Neutral Medium 
(CNM) are in front of this diffuse X-ray emission. This 



amount of absorbing material has to be evaluated in order 
to calculate the strength of the photoelectric absorption. 
We account for the different characteristic scale heights of 
the neutral species by using the velocity integrated maps 
defined in Sect. 13.1. "21 

3.2.3. The extragalactic component 

The extragalactic X-ray background is composed of su- 
perposed em ission of p oint sources, in particular of 
AGNs (e.g. iTozzi et al.l HoOl). Its spectrum can be 
approximated by a power law {E~^). We assume 
an avera ge sp e ctral index of F = 1.5 according to 
Hasingcr et al.l tOO'^. Its ROS A T C- band intensity 
was determined by Barber et al.l l(l99(j ): they derived 
/extra = 228 • 10~^ cts arcmin"^ with an uncertainty 
of 90- 10~^ cts arcmin^^. Throughout the fit procedure 
(see Sect. I3.3|l we will fix /extra to both values normal- 
ized to the ROSAT C-band and their corresponding values 
in the other ROSAT energy bands respectively. We assume 
that all neutral interstellar matter along the line of sight 
measured in the Leiden/Dwingeloo Hi survey (iVn^c) is 
in front of the extragalactic sources. Also this amount of 
photoelectric absorbing material is constant in the evalu- 
ation of the soft X-ray radiative transfer calculations. 

This leads to the radiative transfer equation for soft 
X-rays introduced bv lKerp et al.l 



LHB 



/ 



HALO • e 



/extra ■ e 



(2) 



Here, (Th denotes the effective photoelectric absorp- 
tion cross section of the Galactic halo, while CTo repre- 
sents that for the extragalactic background. While we 
can safely assume that all Hi observed towards a distant 
AGN is enclosed within the velocity range — 450kms^^ < 
''^LSR < 450 km s""'^, we have to question this assumption 
with respect to the Milky Way. Here we have to iden- 
tify the best fit velocity range by comparing the Hi and 
the diffuse X-ray data quantitatively. To differentiate be- 
tween both absorbing column densities we use -/Vhi.g and 
A^Hi,h from here on. Equation |21 has four free parame- 
ters (/lhb, /halo, /lhb, /halo) and two fixed parameters 
(/extra, r). The temperatures of the LHB (Tlhb) and of 
the Galactic halo (Thalo) are connected to Eq. |21by the 
source term S{E) in Eq.Q 7Vhi,o and iVn i^h are fixed pa- 
rameters as well because by optimizing the free parameters 
both A^Hi values are fixed. To determine the four free pa- 
rameters, we use the information of the ROSAT R1-, R2-, 
C- and M-Band in simultaneous fit procedures of X-ray 
energy band ratios and intensities vs. A^Hi,h- 

3.3. Modeling 

In order to evaluate the radiative transport equation (Eq. 
121) for all ROSAT energy bands and band ratios simulta- 
neously, we developed a two-step fit procedure. The first 
step (one dimensional approach) gives us the opportunity 
to estimate the start-values for the fit parameters. The 
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second step {two dimensional approach) is used to derive 
the spatial distribution of the X-ray intensities correctly. 
Moreover, we can optimize the fit parameters within the 
framework of the two dimensional fit procedure. 

3.3.1. The one dimensional approach 

The radiative transport equation (see Eq. gives us the 
possibility to derive fit parameters, namely the intensities 
of the two free emission components. In practice we calcu- 
late so-called scatter diagrams (see Fig.^for the C-band 
scatter diagram) for the ROSAT C- and M-band and the 
energy band ratios R1/R2 and C/M as functions of A'hi- 
The energy band ratios are a sensitive measure for the 
temperature of both X-ray emitting plasmas. 



1 




h . J! 



Fig. 1. Scatter plot for the ROSAT C-band. At high col- 
umn densities (> 5 • 10^" cm^^) the X-ray intensity shows 
an asymptotic behavior. This represents only the fore- 
ground emission of the LHB. The solid line shows the 
best fit model which satisfies all ROSAT energy bands and 
band ratios simultanously (see Fig. [3 El El as well). Other 
models (i.e. with varying halo- and/or LHB-temperatures) 
are not distinguishable from each other. 



The first free parameter value we derive is the C-band 
intensity of the Local Hot Bubble (/lhb ) ■ The data points 
in the C-band scatter diagram (Fig.^l show an asymptotic 
behavior at high column densities (> 5 • 10^°cm~^). Here 
we observe only the X-ray foreground emission because 
the entire background emission of /halo and /extra is 
absorbed by the high column density in front of both 
components; therefore, the first estimate is /lhb = 450 • 
10~^ cts s~^ arcmin"^ which is consistent with the values 
given bv lSnowden et alJ l)l998|) . 

The /halo is determined at the low column density 
part of the C-band and full M-band scatter diagram. In 
the ROSAT C-band we observe the superposed emission of 
/lhb and /halo while in the ROSAT M-band (see Fig.^J 
we have to include the diffuse X-ray emission attributed 
to the extragalactic background. 
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Fig. 2. Scatter plot for the ROSAT M-band. The dotted 
lines correspond to models where the halo plasma temper- 
ature is 10^-^ K (upper curve) and 10^-^ K (lower curve). 
The dashed lines correspond to models where the LHB 
temperature is 10^'^ K (upper curve) and 10^'^ K (lower 
curve). As in Fig.^the solid line represents the best fit 
model. 
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Fig. 3. ROSAT C/M-band ratio vs. iVni. The dotted hues 
correspond to models where the halo plasma temperature 
is 10^ '^ K (upper curve) and 10^'^ K (lower curve). As in 
Fig.n the solid line represents the best fit model. Note, 
that the dashed lines, representing the variation of Tlhb 
in Fig. El are not suitable to evaluate Tlhb in the C/M- 
band ratio. 



Considering only the C-band scatter diagram (Fig.^, 
an ambiguity concerning the plasma temperatures re- 
mains. Variations of the plasma temperature Thalo yield 
insignificant variations in the fitted curves. To constrain 
Tlhb and Thalo we plot the R1/R2 and C/B band ratio 
versus iVnLh (see Fig.|31&EJ. Regarding these, we see large 
variations in the C/M Band ratio vs. -/VHi,h plot, when 
we vary the halo plasma temperature in the same way as 
for the C-band scatter diagram (see Fig.^. Variations of 
Tlhb cause significant variations in the R1/R2 band ra- 
tio vs. A^Hi.h (Fig.Q. An immediate impact on the energy 
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Fig. 4. ROSAT Rl/R2-band ratio vs. Nm- The dashed 
hnes correspond to models where the LHB temperature 
is 10^'^ K (upper curve) and 10^'^ K (lower curve). As in 
Fig-d the solid line represents the best fit model. Note, 
that the dotted lines, representing the variation of Thalo 
in Fig. O are not sensitive to evaluate Thalo in the C /M- 
band ratio. 



band ratios vs. A^Hi.h is seen if we vary /halo a-nd /lhb 
individually. Because a variation of only one free param- 
eter has an immediate effect, all four scatter diagrams, 
we have a sensitive tool to find a consistent solution for 
all four parameters simultaneously. To illustrate this pro- 
posal we plot a set of temperature models for Thalo in 
each of the four scatter diagrams while leaving the other 
parameters /lhb , /halo and Tlhb fixed to the best fit pa- 
rameters given in Sec. 13.4.21 Since the C-band and Rl /R2 
band ratio scatter diagrams are not sensitive to tempera- 
ture variations of the parameter Thalo the M-band and 
C/M band ratio scatter diagrams show up with significant 
deviations from the best fit situation marked by the solid 
line. 

In Fig.^ we demonstrate the variation of the Tlhb 
on the R1/R2 band ratio scatter diagram. While the C- 
and M-band scatter diagrams as well as the C/M ratio 
scatter diagram are not sensitive measures for Tlhb we 
find that the R1/R2 band ratio scatter diagram allows 
one to determine Tlhb very accurately. 

The one-dimensional approach performed here indi- 
cates that the field-averaged values of the four free param- 
eters can be accurately determined simultaneously. The 
Rl /R2 band ratio scatter diagram is a measure for Tlhb ■ 
The C/M ratio scatter diagram is a measure for Thalo- 
The C-band scatter diagram determines /lhb and /halo j 
while the M-band scatter diagram is a measure for /halo 
and Thalo- 

These four parameters form the smallest set of param- 
eters which are necessary to model the X-ray intensity 
distribution in the R1-, R2-, C- and M-band. The one 
dimensional approach leads to the following initial fit pa- 
rameters: 



/lhb = (450±90) • lO^^ctss^^arcmin"^ 
/halo = (1200± 120) - lO^^ctss^^arcmin-^ 
/extra = (230 ± 90) • 10"*^ cts s"^ arcmin"^ 



Thalo = lO^^^^iK. 

3.4. The two dimensional approach 

Extending our analysis of the diffuse X-ray intensity dis- 
tribution from one to two dimensions allows us to disclose 
systematic sources of uncertainties hidden by the scatter 
of the X-ray counts. Namely variations of /lhb and /halo 
as functions of I and b on small as well as on large angular 
scales. 

In our two dimensional approach we evaluate Eq.|21in a 
pixel-by-pixel manner. For simplicity we assume that the 
four parameters are constant across the whole sky region 
we study. 

Equipped with the initial parameters, we calculate 
model maps for the C- and M-band and compare them 
to the observed X-ray distribution in these bands. This is 
done by subtracting the model maps from the observed 
ones and dividing the difference by the ROSAT uncer- 
tainty maps. The result is a deviation map which shows the 
discrepancy of observation and model in units of the sta- 
tistical standard deviation. These deviation maps (Fig.O 
reveal that most of the observed ROSAT C- and M-band 
data are fully consistent with the corresponding modeled 
maps. While the ROSAT M-band data matches with the 
modeled M-band (within the uncertainties), the C-band 
shows some large scale deviations between the observed 
and the modeled X-ray distribution. Especially, the intra- 
group gas of the M81/M82 group of galaxies leads to an 
expected deviation in the ROSAT C-band (see Fig.O up- 
per panel) . This is caused by the admixed H i line emis- 
sion of the intra-group gas and the Galactic hydrogen. 
The supernumerary Hi column density is seen as an X- 
ray shadow in the modeled ROSAT C-band located at the 
correct position on the sky. In this sense, we can test the 
two dimensional approach spatially. 

Towards the low- and high-latitude regions of the maps 
the observed and the modeled maps agree quantitatively 
(Fig.lSJ. Here, the dynamic range in X-ray absorbing Hi 
column density reaches extreme values. The quantita- 
tive agreement of these areas indicates that the general 
approach is successful. However, the differences between 
both map s are pronounced towards the HVC-complex C 
ijKerp et al . 1999) and the Lockman Hole. Here, we focus 
on the Lockman Hole. 

3.4.1. Testing the two dimensional correlation 

To optimize the initial fit parameters (see Sec. I3.3|l we 
compute histograms of the deviation maps and obtain the 
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Fig. 5. upper panel: The color image shows the observed 
ROSAT C-band while the contours correspond to the de- 
viation map. The solid contours encircle regions where the 
modeled intensities are too faint, the dotted contours en- 
close regions where the model is too bright. Intra-group 
gas of the M81/M82 group of galaxies within the spec- 
tra of the Leiden/Dwingeloo Hi survey produces the con- 
tour line around {l,b) « (140°, 40°). HVC complex C is 
located at the bottom, the Lockman hole is surrounded 
by the dotted contour line on the right hand side, lower 
panel: The color image shows the modeled ROSAT M- 
band. The superimposed contours of the corresponding 
deviation map are almost nonexistent, indicating a quan- 
titatively good model within the uncertainties. Intensities 
are in 10~^ctss~^ arcmin"^. 



mean (/i) and variance (ct). In case of a perfect quantita- 
tive correlation between the observation and the model we 
expect a mean of ^ = and a variance of ct = 1 . This can 
be achieved by small variations of the free fit parameters 
and slight changes in the H i column density of the halo 
component (i.e. the velocity interval of the Hi data). It 
turns out that the quality of the histograms increases when 
the intensity of the extragalactic background emission is 
tuned to (170 d z 90) ■ 10~^ ct s s~^ a rcmin^^ which is fully 
consistent with lBarber et all l)l996|) within the uncertain- 
ties. After seven iteration steps no further improvement 
can be achieved. The statistical parameters in the C- and 



M-band after the optimization are fj,c — 0.13, ctc — 1.56 
and /iM — —0.58, au = 1.12 respectively. The final C- 
band model map is shown in Fig.|Hl for comparison with 
the observed C-band intensity distribution (Fig.El). 

3.4.2. The role of Nui 

In order to test which H i components are responsible for 
the major absorption we use different A^hi maps (see Sect. 
12. 1|) . Obviously the smaller the velocity integration in- 
terval the lower the A^hi column density. However, the 
^lsrI < 25kms~^ map does not incl ude the IVC shad - 
ows - most obvious for IVC 135-^52-45 l(Wei6 et al.lll999t) . 
The X-ray excess emission close to HVC complex C is still 
present, as well as the deficit of attenuating H i towards the 
Lockman Hole. Figures IHIIZI and |H1 show three model maps 
of the ROSAT C-band intensity distribution. The maps 
differ in the amount of H i in the halo component. We cal- 
culate the maps for three different LSR velocity intervals: 



-25kms ^ < ULSR < -l-25kms \ 



-50 kms 



< t^LSR < 

-1-50 km s"^ and -70kms"^ < wlsr < -h70kms"\ The 
model map to the smallest velocity interval has a mean 
intensity of Jc = (960 ± 360) • 10"^ ctss^^ arcmin"^ ^pj^g 
model map to the medium velocity interval differs by 11 % 
and the model map to the largest interval by 18% (see 
Fig. El - IHll. From these values we conclude, that most 
of the CNM and WNM is already included within the 
smallest velocity range, which is accordingly responsible 
for the dominant fraction of the absorption. Moreover, the 
CNM is not present at the positions of low column density 
l|Lockman et al.iri986(l and therefore the WNM dominates 
the absorption of the X-ray radiation. This statistical ap- 
proach leads to the optimized fit parameters: 



^HALO 
^EXTRA 



(350 ± 90) • 10"*^ cts s^i arcmin"^ 
(1380 ± 120) • 10"^ cts s"^ arcmin-2 
(170 ± 90) • 10^^ cts s"^ arcmin-2 



Tlhb = 106°±«iK 
Thalo = 106-2±oiK. 



In comparison to iKuntz fc SnowdenI l|200Cl() we only 
need one single X-ray emitting component for fitting the 
Galactic halo emission. The temperature wc derive is in- 
between t he temperature of the soft and hard component 
derived bv lKuntz fc SnowdenI l|200Cll) . The temperature for 
the LHB component is slightly lower in our case. 

It turns out that the entire region around the Lockman 
Hole {l,b) = (150°, 53°) is correctly modeled down to the 
data uncertainties. We can account for the X-ray intensity 
distribution on a large sky region (60° x 60°) whereas the 
Lockman Hole itself is too bright in the model. The ques- 
tion why the model first fails at this prominent position 
will be discussed in the forthcoming section. 

One of the main results is that Tlhb < Thalo- This 
finding is obvious in the R1/R2 band ratio vs. A^Hi,h- 
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Fig. 6. Modeled ROSAT C-band map. The amount of H i 
corresponds to the velocity interval of |wlsr| < 25kms^^ 
which does not include the IVC regime. Therefore, 
the X-ray shadow produced by IVC 135+52-45 (com- 
pare to Fig.[Sl upper panel) is absent. Intensities are in 
10~^ cts s~^ arcmin"^. 




Fig. 7. Modeled ROSAT C-band map for the H i veloc- 
ity interval of |ulsr| < 50kms~^. The mean intensity is 
about 11% smaller than the one of Fig.|Hl Intensities are 

in 10^® cts s^^ arcmin^^. 

Because of the strong dependence of the photoelectric ab- 
sorption on energy (cr oc E~^), the soft X-ray photons are 
much stronger attenuated than the hard ones, leading to 
an apparent hardening of the X-ray spectrum proportional 
to A^Hi- This implies that the higher A^Hi,h the lower the 
R1/R2 band ratio of the /halo source component. Figure 
0|shows the contrary behavior. The higher -/VHi,h the larger 
the R1/R2 band ratio. This finding can only be explained 
by a cooler LHB plasma with respect to the Galactic halo 
gas. Using the best fit plasma temperatures we derive 
an important contribution of the Galactic halo plasma 
to the observed M-band intensities. The /halo is about 
4 times brighter than /lhb within the ROSAT C-band, 
while /halo — |/extra in the ROSAT M-band. /halo 
is the dominant source in the ROSAT C- and M-band. 
We assumed the intensity of the Galactic halo plasma to 




Fig. 8. Modeled ROSAT C-band map for the Hi veloc- 
ity interval of |wlsr| < 70kms^^. The mean intensity is 
about 18% smaller than the one of Fig.|Sl Intensities are 
in 10~^ ctss~^ arcmin^^. 

be constant across the field of interest which is consistent 
within the uncertainti es. The derived fluctuations in halo 
intensity suggested bv lSnowden et alJ l(l99 4a) may be due 
to the Hi velocity intervals, which differ from our choice. 
These fluctuations do not appear in our model, since we 
choose the H i velocity intervals as described above. 

4. Specifying the ISIVI conditions in the Lockman 
Hole 

4.1. Molecular or ionized hydrogen? 

As we pointed out in the prior sections, most of the sky 
region investigated can be explained by our model. The 
deviating portions of the sky decompose into two types: 
regions modeled too faint and those modeled too bright. 
Here we concentrate on the Lockman Hole, which belongs 
to the latter case (see Fig.|S| upper panel). This is due 
to a lack of absorbing material. The additional hydrogen 
column density required to remove the deviation can be 
estimated to N^'^'^ « (5 ± 3) • 10^^ cm'^. 

Since most of the sky region is properly fit by our 
model, we argue that the ISM towards the Lockman Hole 
has a different composition to the surroundings. Either 
the ISM at this particular position contains a molecular 
gas phase or an ionized gas phase, both of which are de- 
tectable as X-ray shadows but not quantitatively traced 
by the Hi data. The Hi column density inside the bor- 
dering contour lines in Fig.|Sl (upper panel) is less than 
lO^^cni^^. If we consider a molecular gas phase, a mini- 
mum H I column density is needed to shield this molecu- 
lar phase. Herbstmeier et al.l l|l993() . for example, claimed 
that an Hi column density of at least 3.7 • 10^" cm^^ 
is needed to allow the formation of molecular hydrogen. 
From this, we conclude that a molecular gas phase is un- 
likely to exist. The remaining option is an ionized state of 
the hydrogen towards the Lockman Hole. Previous studies 
of the interstellar H„ emission towards the Lockman Hole 
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ijHausen et al.ll2002|) yield iVHn ~ 2 • 
the same order of magnitude as N^^'^. We perform a cross- 
check of the ionization structure towards the Lockman 
Hole via far-ultraviolet absorption line measurements by 
FUSE. 

4.2. Far-Ultraviolet absorption line data 

In order to investigate the physical conditions of the inter- 
stellar medium in the direction of the Lockman Hole we 
make use of Far-Ultraviolet (FUV) absorption line data. 
The FUV range bluewards of the Lyman a line and red- 
wards of the Lyman break provides a wealth of spectral 
diag nostics to study inte r stellar gas in its vario us phases 
(e-g.- ISavage et alJ l)2003l) : lRichter et all l)2001a|) ). A large 
number of FUV absorption line spectra for almost every 
direction in the sky is publicly available from the MAST 
data archive of the Far Ultraviolet Spectroscopic Explorer 
(http : / / archive . stsci . edu/f use/). For the purpose of 
studying the Lockman Hole field extragalactic background 
sources such as QSOs and AGNs are favorable in order to 
ensure a correct continuum placement for the various ab- 
sorption features, and to include Milky Way halo gas that 
may contain a significant fraction of the X-ray absorbing 
neutral and ionized gas in this general direction. 

4.3. NGC 3310 

From the extragalactic sight lines in the general direction 
I ~ 150°, b ~ 50° for which good S/N FUSE data is 
publicly available, we have chosen the one closest to the 
Lockman Hole: NGC 3310 [l = 156.6°, b = 54.1°). The raw 
data was extracted using the standard FUSE data reduc- 
tion pipeline (CALFUSE v. 2. 1.6), and was then further 
reduced in a way similar to the data processing procedure 
described in Wakkcr ct al. (2003). The velocity calibra- 
tion was made by comparing absorption lines from vari- 
ous atomic and molecular lines with existing H i 21 cm-line 
data (see below) . The FUSE data have an average signal- 
to-noise ratio (S/N) of ~ 14 per resolution element, and a 
spectral resolution of 20 kms"^ (FWHM). 

4.3.1. Neutral and weakly ionized gas 

In Fig.ini (uppermost panel) we have plotted the 
Green Bank H i 21 c m-line spectrum towards NGC 3310 
llMurohv et alJ ll99(j^ on a LSR velocity scale. The 21 cm- 
line spectrum shows at least four distinct velocity com- 
ponents: (1) a broad local Galactic component centered 
at -11 kms-i with log N(Hi) = 19.63 ± 0.02, (2) a very 
narrow local Galactic component at —17 kms~^ having 
log N(Hi) = 19.54 ± 0.01, (3) a relatively strong broad 
component located at —47 kms""'^ with log N(Hi) — 
19.71 ±0.01, re l ated t o the Intermediate- Velocity Arch (IV 
Arch; IWakkeJ l|200J) ) in the lower Milky Way halo, and 
(4) a weaker, second IV Arch component at —65 kms~^ 
with log N(Hi) = 19.11 ± 0.01. In the lower four pan- 
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Fig. 9. From top to bottom, the panels show the Hi 
21 cm-line spectrum towards NGC 3310 as measured by 
the Green Bank Telescope, and the absorption profile of 
Oi, An, Pii, and Ha P(2),4-0 as measured by FUSE. 
The dotted lines indicate the LSR velocities —5, —45, and 
-70kms"i (see Sect. EXTl for details). 



els of Fig.|51we have plotted the FUSE absorption profiles 
of Oi A936.6, An A1048.2, Pii A1152.8, and Ha P(2),4-0 
A1053.3 in direction of NGC 3310. 

Only three velocity components (fitted by Gaussian 
absorption profiles) are visible in absorption, centered at 
approximately —5, —45, and —70 kms~^, in Fig. |51 in- 
dicated by dotted lines. Possibly, the two local Galactic 
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Fig. 10. The Ovi A1031.9 absorption profile towards 
NGC3310 traces the highly ionized gas (see Sect. 14.3!^ 
for details). 



components, in the 21 cm-line emission data seen at —11 
and —17 kms~^, are smearing together into one single 
absorption component, given the fact the FUSE data has 
a spectral resolution of only ~ 20 kms~^. However, it 
is also possible that the narrow —17 kms~^ 21 cm-line 
component relates to gas that partly fills the 21.'0 Green 
Bank beam, but is not present along the pencil-beam 
FUSE spectrum towards NGC 3310. Since beam smearing 
effects can introduce enormous uncertainties when com- 
paring 21 cm- line data with FUV absorption line data for 
low-column density clouds, it is important to look for pos- 
sibilities to study the structure of the neutral interstellar 
gas towards the Lockman Hole based on the FUSE data 
alone. Unfortunately, the Hi Lyman series in the FUSE 
bandpass cannot be used for this purpose, since these lines 
are heavily saturated at this column density level. 

The Ar i A1G48.2 (Fig.lHI third panel) indicates that the 
bulk of the neutral material resides in the IVC component 
near —45 kms~^. And also the H2 absorption (lower-most 
panel) suggests that most of the diffuse molecular gas is 
situated in the IV Ar ch rather than in th e local Galactic 
component (see also iRichter et al.l l)2003l) for details on 
the H2 spectrum). The best available tracer for the neu- 
tral gas in the FUSE bandpass is Oi. The Oi column 
density is tightly bound to the H I column density, since 
both neutral hydrogen and neutral oxygen have the same 
ionization potential (13.6eV), and their abundance is cou- 
pled by charge-exchange reactions. In addition, oxygen is 
not strongly depleted in dust grains. Unfortunately, O i 
absorption towards NGC 3310 is saturated for all avail- 
able lines at A > 930 A, and only the relatively weak 
Oi A936.6 can be used to roughly estimate the Oi col- 
umn density. This line is shown in Fig.|51 second panel. 
As is clearly visible, the O i absorption does not reach the 
zero-flux level in any of the three absorption components, 
which (together with the total width of the composite ab- 
sorption structure) allows upper limits to be placed on 
the individual equivalent widths and column densities by 
modeling the component structure with Gaussian absorp- 
tion components. Such a model also has to account for the 
H2 P(3) absorption that partly blends the local Galactic 



Oi absorption component near -|-50kms"^ (in Fig.|^ in- 
dicated with 'b'). We have done such modeling for the O i 
A936.6 line, and find that the O i absorption in the local 
Galactic component at —5 kms~^ must be < 90 mA {3a 
upper limit), while for the two IVC components at —45, 
and —70 kms~^ we find equivalent width limits of 115 and 
85mA. Assuming that b = lOkms"^ (or greater), these 
equivalent width limits correspond to logarithmic column 
density limits of < 16.02, < 16.57, and < 15.94 for Ol, 
and < 19.28, < 19.83, and < 19.20 for Hi, respectivel y, 
if [0/H] is solar (log (0/H)q = -3.26 llHolwegeill200ll) 1. 
Thus, the FUSE data suggest that the Hi column den- 
sity in the local Galactic gas towards NGC 3310 is by a 
factor of ^ 2 lower than indicated by the beam-smeared 
21 cm-line data pointed in this direction, while the 21 cm- 
line matches the FUSE data for the two IVC components. 
However, unresolved sub-components (particular in the lo- 
cal Galactic component) may be present, in which case 
these limits could be incorrect. 



What does the FUSE spectrum of NGC 3310 teU us 
about the ionization conditions in the interstellar gas 
in the direction of the Lockr nan Hole? As has been 
demonstrated in earlier studies l|Richter et al ] l2001b^ . the 
Pii/Oi column density ratio is a good indicator for the 
ionization structure in interstellar clouds. P 11 has an ion- 
ization potential of 19.7 eV, and lives in regions where O 
and H are both neutral and ionized. Also for phosphorus, 
dust depletion is not considered to be important. Thus, 
assuming an intrinsic P/0 ratio, the Oi and Pii column 
densities can be used to determine lower limits for the ion- 
ization fraction, /ion =H+/(H°-|-H"'"), in the weakly ion- 
ized gas phase (for ionization energies < 19.7 eV). The 
weak P II A1152.8 absorption (Fig.lHl fourth panel), can be 
reproduced by a three-component Gaussian fit, for which 
the three velocity components have been fixed at —5, —45, 
and —70 kms^^. Using this method, we derive equivalent 
widths of 18 ±4 mA for the —5 kms~^ component, 60 ±8 
niA for the —55 kms^^ component, and 11±4 mA for the 
—70 kms^^ component. Assuming a b value of 10 kms~^, 
these equivalent widths correspond to column densities, 
logN(Pii), of 12.85 ±0.10, 13.51 ±0.09, and 12.63 ±0.14, 
respectively. Assuming that all three velocity components 
hav e solar phosphorus-to-o xvKen abundance ratios (see 
also lRichter et aObflOlalbh . the data for O and P yield 
/ion > 0.22 and N(H)==N(H i)+N(H+) ^ 8.66 x lO^^cm'^ 
whereas no useful limits can be obtained for the local 
Galactic component and the second IVC component, since 
their P/0 column density ratios would allow every degree 
of ionization (i.e., /ion > 0.0). The FUSE data therefore 
indicates that approximately a fifth of the total hydrogen 
column density in the main IVC component is in ionized 
form; very likely, the ionization in the local Galactic com- 
ponent is much higher, but the limitations for the FUSE 
O I data unfortunately does not allow the derivation of a 
quantitative limit. 
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4.3.2. Highly ionized gas 

The FUSE data also yields information about the highly 
ionized gas phase for ionization energies > 100 eV through 
the two resonance lines of five-times ionized oxygen, O vi, 
situated at 1031.9 and 1037.6 A. Ovi traces collisionally 
ionized, hot gas at temperatures T ~ 2 x 10^ K (thus lower 
than the X-ray emitting gas), or gas that is photoionized 
by a hard UV background (or a combination of both) . In 
Fig-El we have plotted the Ovi A1031.9 absorption pro- 
file towards NGC3310 (the weaker Ovi A1037.6 cannot 
be used due to blending efi'ects). The main analysis of the 
Ovi A1031.9 abs orption toward s NGC 331 has already 
been p resented in lWakker et al.l l)2003|) and I Savage et al.1 
l|2003(l . Ovi absorption is smeared over a broad range 
from —200 to -1-90 kms^^, thus spanning a velocity range 
that covers all three velocity components seen in 21 cm- 
line emission. The total integrated O vi column density is 
log N(0 vi) = 14.56 ±0.02. This column density is slightly 
(but not significantly) higher than the Ovi column den- 
sities in other dire ctions outside the Lockman Hole (see 
I Savage et alJ l2003l) . This is not too surprising though, 
since the pathlength through the Lockman Hole is too 
short to contribute significantly to the total O vi column 
density at the expected low volume densities. Assuming 
that the fraction of oxygen residing in the O vi phase is 
< 0.2, the total gas column density of this highly ionized 
gas pha se is < 3.3 x 1 0^^ cm^^ for solar oxygen abun- 
dances l)Ho1wegeT< 12001 . The highly ionized gas compo- 
nent traced by Ovi therefore cannot contribute signifi- 
cantly to the total hydrogen gas column density in direc- 
tion of the Lockman Hole. 

5. Conclusion 

The analysis performed shows that the main absorption of 
soft X-rays is due to the WNM. It turns out that a single- 
temperature plasma in the Milky Way halo is sufficient 
to reproduce the X-ray distribution on the Lockman Hole 
field (60° X 60°) very well. Hence, a patchy Milky Way 
halo is not likely. Furthermore, the applied fit procedure 
(simultaneously for four ROSAT energy bands) is very 
sensitive to temperature. We showed that the temperature 
of the LHB plasma is lower than the temperature of the 
halo plasma. The X-ray/Hi-analysis of the Lockman Hole 
itself gives evidence that a large fraction of the attenuating 
ISM may reside in form of ionized hydrogen. 

While the FUSE data of NGC331G impressively 
demonstrates the complex multiphase structure of the gas 
in the general direction of the Lockman Hole, a precise 
quantitative analysis of the gas properties turns out to be 
unfeasible due to the complexity of the velocity structure 
of the gas that cannot be resolved with the FUSE instru- 
ment. The data shows that the majority of the neutral 
and molecular gas is situated in the IVC component in the 
lower Milky Way halo rather than in the disk component. 
Due to the data limitations, we are not able to obtain pre- 
cise ionization fractions, but the FUSE data of NGC 3310 



are at least consistent with the idea of ionized hydrogen 
in the Lockman Hole. While the ionization fraction is not 
well constrained by the FUSE data the total amount of 
ionized gas is sufficiently high to overcome the apparent 
discrepancy between the Hi and X-ray observations. 
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